Dynamic compartmentalization of bacteria: accurate division in E. coli. 
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Positioning of the midcell division plane within the bacterium E. coli is controlled by the min 
system of proteins: MinC, MinD and MinE. These proteins coherently oscillate from end to end of 
the bacterium. We present a reaction-diffusion model describing the diffusion of min proteins along 
the bacterium and their transfer between the cytoplasmic membrane and cytoplasm. Our model 
spontaneously generates protein oscillations in good agreement with experiments. We explore the 
oscillation stability, frequency and wavelength as a function of protein concentration and bacterial 
length. 



The subcellular spatial and temporal organization of 
bacterial proteins is largely unknown. Already, the spa- 
tial distribution of proteins on the cytoplasmic membrane 
of bacteria are known to be important for chemotaxis 
and for DNA replication . Improving our understand- 
ing of how this supra-molecular organization of proteins 
affects bacterial function represents a considerable exper- 
imental and theoretical challenge. In contrast to nucle- 
ated eukaryotic cells, no large organelles are present in 
the bacterial interior (cytoplasm) and no active trans- 
port mechanisms such as molecular motors are known to 
function there. However, recent video microscopy of fiuo- 
rescently labeled proteins involved in the regulation of E. 
coli division have uncovered coherent and stable spatial 
and temporal oscillations in three proteins: MinC, MinD, 
and MinE §, g |, |, |, §. The proteins oscillate from 
end to end of the bacterium, and move between the cy- 
toplasmic membrane and the cytoplasm. These min pro- 
teins select the site for the next bacterial division || [l0| . 
Despite a wealth of phenomenological detail, no quanti- 
tative models have been developed of how the min pro- 
teins organize into oscillating structures. Understanding 
the self-organized patterns involved in bacterial division 
processes can give us insight into how a bacterium can 
dynamically compartmentalize itself. 

We focus on E. coli, a commonly-studied rod shaped 
bacterium, approximately 2 — 6 fim in length and around 
1 — 1.5/im in diameter. Each E. coli divides roughly ev- 
ery hour, depending on the conditions — first replicating 
its DNA then dividing in half to form two viable daugh- 
ter cells. The MinCDE oscillations are known to per- 
sist even when protein synthesis is suppressed [||, and 
DNA replication and septation occur even without the 
min proteins. Hence the min system can be studied inde- 
pendently of the other division processes. Efficient divi- 
sion requires many processes, including DNA replication, 
MinCDE oscillations, and the actual septation process. 
Septation initiates with a contractile polymeric "Z-ring" 
of a tubulin-homologue FtsZ that forms just underneath 
the cytoplasmic membrane. The FtsZ septation rings 



largely avoid guillotining the DNA-containing nucleoids 
independently of the min system . This "nucleoid oc- 
clusion" serves as a complementary control mechanism 
for accurate cell division. The role of the min system ap- 
pears to be to restrict the Z-ring to midcell. This reduces 
the production of inviable nucleoid-free minicells which 
occur when the cell divides too close to the cell poles. If 
the min system is genetically knocked out, 40% of divi- 
sions lead to inviable minicells || — a sizeable drain on 
bacterial resources. 

The study of deletion mutants has made the phe- 
nomenological roles of the individual min proteins clear. 
MinC associated to the cytoplasmic membrane locally 
inhibits assembly of the contractile Z-ring, but remains 
cytoplasmic and largely inactive in the absence of MinD 
H. MinD binds MinC and recruits it to the cytoplas- 
mic membrane [l2|. MinE drives MinD away from 
the bacterial midplane, and hence allows a contractile 
ring to form only there. Without MinE, the membrane- 
bound MinC/MinD block Z-ring formation everywhere, 
inhibiting division, and resulting in the formation of long 
filamentous cells §|, §. Without MinC, Z-ring formation 
cannot be inhibited anywhere and inviable minicells are 
produced. Without MinD, neither MinC nor MinE are 
recruited to the cytoplasmic membrane and so have re- 
duced effect. 

With normal levels of MinC, D, and E, a remarkable 
oscillatory dynamics is seen ||, |J, Hi |f| @> Hi- First the 
MinC /MinD accumulate at one end of the bacterium on 
the cytoplasmic membrane. Then MinE forms a band at 
midcell which sweeps towards the cell pole occupied by 
the MinC/MinD, ejecting the MinC/MinD into the cyto- 
plasm as it goes. The ejected MinC/MinD then rebinds 
at the other end of the bacterium. When the MinE band 
reaches the cell pole, it disassociates and reforms at mid- 
cell. The entire process then repeats towards the opposite 
cell pole. The oscillation period is approximately 1 — 2 
minutes, so many oscillations occur between each bacte- 
rial division. The dynamics minimizes the MinC/MinD 
concentration at midcell, thereby allowing the Z-ring and 
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the subsequent division septum to form there. 

Until recently || , qualitative models of the min system 
involved unidentified midcell topological markers (see, 
e.g., Q). This letter puts forward the first quantitative 
self-organized model that describes much of the intri- 
cate phenomenology of accurate division site placement 
in E. coli, and does so using only the diffusive motion and 
interactions of the min proteins. The essence of our ap- 
proach is to describe the MinCDE dynamics by a set of 
coupled reaction-diffusion equations. Experimental re- 
sults indicate that the oscillatory protein dynamics is 
unaffected if new protein synthesis is blocked || . Accord- 
ingly we employ a model that conserves the total number 
of each protein type. Strikingly, this model possesses a 
linear Turing-like (Hopf) instability |l5) despite the 
absence of mechanisms such as internal reactant produc- 
tion or external feed that have normally been required 
to model Turing patterns [Of course energy input 

in the form of ATP is required to sustain the oscillations 
within a bacterium.] As we will see, the resulting pro- 
tein oscillations mark the midcell with a minimum of the 
time-averaged concentration of MinC /MinD and with a 
corresponding maximum of MinE. 

Our starting point is a set of four coupled reaction- 
diffusion equations describing, respectively, the densities 
of MinD on the cytoplasmic membrane (fid), MinD in 
the cytoplasm (pd), MinE on the cytoplasmic membrane 
(p e ), and MinE in the cytoplasm (pe)- 
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Following the observation in Refs. j|, |[ that the MinC 
dynamics simply follows that of the MinD, we do not 
model the MinC field explicitly. We consider the varia- 
tion of density along the long bacterial axis, tracking the 
local rates of change of the densities stemming from dif- 
fusion and from transfer between the cytoplasmic mem- 
brane and the cytoplasm. Zero flux "closed" boundary 
conditions are imposed at both ends of the bacterium. 
The total amount of MinD and MinE, obtained by in- 
tegrating pd + pD and p e + pe over the length of the 
bacterium, is explicitly conserved by our dynamics. 

By reducing the min protein dynamics to a set of de- 
terministic Id rate equations we neglect fluctuation ef- 
fects. Given that the number of min molecules in each 
cell is rather small (around 3000 for MinD jl7j and 170 for 
MinE Q) these fluctuations could be important. While 
some fluctuation effects are evident experimentally, such 
as an occasional mid-cycle reversal of the direction of 



MinE band propagation g , on the whole bacterial oscil- 
lations appear to be amazingly regular [^J. Our contin- 
uum coarse-grained approach captures the essence of the 
protein dynamics and explains the self-organized aspects 
of the MinCDE oscillations. 

In the first reaction terms in Eqs. (||) and o\ de- 
scribes the spontaneous association of MinD to the cyto- 
plasmic membrane H . MinD is required to recruit MinE 
to the cytoplasmic membrane, but it is an open ques- 
tion whether it is cytoplasmic MinD or membrane-bound 
MinD that is primarily active. A cytoplasmic interaction 
between MinD and MinE has been observed in Ref . [^2) , 
and we are currently only able to obtain the MinCDE 
oscillations by allowing cytoplasmic MinD to recruit cy- 
toplasmic MinE to the membrane, via 03 in Eqs. (|^) and 
(Q). Once on the membrane, MinE drives MinD into 
the cytoplasm. We represent this with 0% in the sec- 
ond reaction terms in Eqs. (0) and (|). Finally, MinE 
will spontaneously disassociate from the membrane, cor- 
responding to (74 in the second reaction terms in Eqs. (|^) 
and (Q). We have not included spontaneous MinD disas- 
sociation or spontaneous MinE association terms, since 
experimentally MinE dominates the MinD disassociation 
and MinD dominates the MinE association. 

Many other reaction terms are possible, however we 
only include the simplest possible "renormalizations" of 
the basic recruitment and release terms, a[ and tr 4 . Ef- 
fectively a[ corresponds to membrane-bound MinE sup- 
pressing the recruitment of MinD from the cytoplasm, 
and 04 to cytoplasmic MinD suppressing the release of 
membrane-bound MinE. We have also set the diffusion 
constants for the membrane-bound MinD and MinE to 
zero. Our results are not qualitatively changed by us- 
ing nonzero values, provided the membrane diffusion con- 
stants remain much less than their bulk counterparts. 

For our simulations we discretized space and time with 
spacings of dx = 8 x 10~ 3 ^m and dt = 1 x 10~ 5 s. We 
have checked that our results are unchanged with smaller 
dx and dt. Densities are measured in molecules per mi- 
cron, and unless otherwise stated we use average densities 
of 1500 pm- 1 for MinD Q and Sbprn" 1 for MinE |JJ. 
The numerical values of our other parameters have not 
been experimentally determined for the min proteins. We 
choose cytoplasmic diffusion constants slightly less than 
the value 2.5 /jm 2 s _1 directly measured for a maltose 
binding protein p9| within the E. coli cytoplasm. Un- 
less otherwise mentioned we use a length of 2pm and 
the following values for the parameters in Eqs. (Q)-^): 
D D = 0.28 pm 2 /s, D E = 0.6 pm 2 /s, a x = 20 s -1 , 
<j[ = 0.028 pm, o- 2 = 0.0063 pm/s, a 3 = 0.04 pm/s, 
(T4 = 0.8 s _1 , and er 4 = 0.027 pm. 

We have analyzed the linear stability of Eqs. (|l|)-(Q) 
p"5| . Testing solutions of the form e xt + l i x with the above 
parameter values, we find a complex \{q) with a posi- 
tive real part that is maximized for q w 1.5/zm , where 
A max = 0.010 ± 0.043i. This indicates the presence of 
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a maximally linearly unstable oscillating mode with a 
wavelength of 4.2 /im and period of 145 s. This finding 
is confirmed by a direct numerical stability analysis of our 
model (not shown). The physical origin of this instability 
lies in the disparity between the membrane and cytoplas- >< 
mic diffusion rates, and also in the slower rate at which £ 
MinE disassociates from the membrane. This ensures A 
that the MinE dynamics lags that of the MinD, setting -i- 
up the oscillating patterns. The existence of the linear v 
instability in Eqs. (|l])-(Q) is crucial, since it means that 
the oscillating pattern will spontaneously generate itself 
from a variety of initial conditions — including nearly 
homogeneous ones. In our simulations we used random 
initial conditions, although identical patterns were also 
observed with asymmetric initial distributions of MinD 
and MinE. The eventual oscillating state is stabilized by 
the nonlinearities in Eqs. (|l|)-(§). At the midcell, this 
oscillating pattern has a minimum of the time-averaged 
MinD concentration — an essential feature of division 
regulation — and a maximum of the time-averaged MinE 
concentration. 




FIG. 1: Space-time plots of the total MinD (left) and MinE 
(right) densities. The greyscale runs from 0.0 to 2.0 times the 
average density of MinD or MinE, respectively. The MinD 
depletion from midcell and the MinE enhancement at midcell 
are immediately evident. Time increases from top to bottom, 
and the pattern repeats indefinitely as time increases. The 
greyscale reference bar spans 100 sec. The horizontal scale 
spans the bacterial length (2 /im). 

Space-time plots of the MinD and MinE concentra- 
tions for a cell length of 2 fim are shown in Fig. [l| In 
excellent agreement with the experimental results, the 
MinE spontaneously forms a single band at midcell which 
then sweeps towards a cell pole, displacing the MinD, 
which then reforms at the opposite pole. Once the MinE 
band reaches the cell pole it disappears into the cyto- 
plasm, only to reform at midcell where the process re- 
peats but in the other half of the cell. These patterns are 
stable over at least 10 9 iterations (10 4 s) — long enough 
for the min system to regulate cell division throughout 
the division cycle of the cell. In Fig. p|, we plot the time- 
averaged MinD and MinE densities as a function of po- 
sition. MinD shows a pronounced dip in concentration 
close to midcell, which allows for the removal of division 
inhibition at midcell. This is in qualitative agreement 



with the experimental data of Ref. ||. MinE peaks at 
midcell, with a minimum at the cell extremities. 
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FIG. 2: The time-average MinD (left) and MinE (right) den- 
sities, {p(x))/pmax, relative to their respective time-average 
maxima, as a function of position x (in um) along the bac- 
terium. 

We also investigated longer filamentous bacteria and 
found a multiple MinE band structure (not shown) . Mul- 
tiple MinE bands always combined into a single MinE 
band in cell lengths shorter than the natural wavelength 
indicated by linear stability analysis. 

The oscillation period as a function of the average 
MinD concentration is shown in Fig. || (left). We find 
a linear relationship indicated by the best-fit line, where 
the period approximately doubles as the MinD concen- 
tration is quadrupled. A linear relationship has also been 
suggested experimentally Q . The period of oscillation as 
a function of cell length is shown in Fig. [3] (right). Below 
lengths of 1.2 [im the bacterium does not sustain oscil- 
lating patterns. For lengths above this minimum, the 
oscillation patterns are stable and the period increases 
with length — as observed experimentally . The peri- 
ods measured from our numerics for cell lengths of 2 /im 
are around 100 s, in good agreement with experiments, 
where periods from 30 — 120 s have been found || . A sin- 
gle MinE band state is stable over a wide range of lengths 
for a given density of min proteins. This provides strong 
evidence that the min system is capable of regulating ac- 
curate cell division over normally occurring cell lengths as 
the cell grows between division events. At longer lengths 
of around 6 [im we observe long-lived metastable states 
with two MinE bands. These multiple bands can survive 
for a thousand seconds or more before decaying into a 
single band. At still longer lengths the two band state 
appears stable; this occurs around 8.4 fim — twice the 
dominant wavelength given by the linear stability anal- 
ysis. This explains why the characteristic wavelength of 
linear stability analysis is rather longer than a normal 
E. Coli bacterium - if the length scale were smaller then 
multiple MinE bands might occur in bacteria of normal 
lengths and proper division regulation would be inhib- 
ited. 

If the MinD concentration is increased or decreased 
beyond the limits shown in Fig. || (left), then the oscil- 
lation amplitude decays, and a uniform steady-state re- 
sults. The stability is mapped out in Fig. m as a function 
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FIG. 3: Left: Plot of the period of oscillation (in seconds) 
against MinD density (in fim^ 1 ), at fixed average MinE con- 
centration of 85 fim" 1 . The solid line is a linear best fit. 
Right: Plot of oscillation period against cell length, for fixed 
MinD and MinE concentrations. Below bacterial lengths of 
1.2/j,m oscillation is not observed. 



of protein concentration. This is consistent with experi- 
ment, where overexpression of MinD suppresses division 
H . Although varying the MinE concentration does effect 
the region of oscillatory instability (as shown in Fig. ^) , it 
did not have a significant effect on the oscillation period. 
This appears somewhat contrary to the results of Ref . || , 
possibly due to the absence of MinE dimerization in our 
model ifj. 



quire that the membrane-associated diffusion constants 
for MinD and MinE are very much less than their corre- 
sponding values in the cytoplasm; and we have mapped 
out the shape of the oscillation regime as a function of 
average MinD and MinE concentration. 

Experimental characterization of reaction rates and 
diffusion constants do not yet severely constrain our 
model. Accurate experimental measurements of oscil- 
lation periods and wavelengths as a function of concen- 
trations of MinD and MinE will provide a stringent test. 
There is also considerable scope for extending our results. 
In subsequent studies we will explore a bulk 3d system 
with discrete particle dynamics and microscopic interac- 
tions between individual protein molecules. This will al- 
low us to explicitly consider the influence of fluctuations 
due to discrete particles, the role of ATPase activity of 
MinD 0, and the effects of MinE dimerization 
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FIG. 4: Filled circles indicate regions of linear instability vs. 
the density of MinD and MinE, where small inhomogeneities 
grew into a periodically oscillating pattern. Open circles in- 
dicate regions of linear stability where small inhomogeneities 
decay into a uniform and static pattern. 

In conclusion, we have introduced a particle- 
conserving reaction-diffusion model that self-organizes 
to form a key regulatory mechanism for accurate midcell 
division site selection in E. coli. The model qualitatively 
agrees with many of the features found in experiments, 
and, in particular, naturally accounts for the oscillatory 
patterns of the min proteins. Already our model leads us 
to make a number of striking predictions: we require that 
cytoplasmic MinD recruits MinE to the membrane; we re- 
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